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Abstract. For the weight function W M (x) = (1 — ||a;|| 2 ) f ', fi > —1, A > and 
fe M a normalizing constant, a family of mutually orthogonal polynomials on the 
unit ball with respect to the inner product 



</,<?> 



/ f{x)g(x)W fl (x)dx + X f Vf(x) ■ Vg{x)W^(x)dx 
Jm d Jm d 



are constructed in terms of spherical harmonics and a sequence of Sobolev 
orthogonal polynomials of one variable. The latter ones, hence, the orthogonal 
polynomials with respect to (■,■}, can be generated through a recursive formula. 



1. Introduction 

The orthogonal polynomials on the unit ball M d of R d have been well understood, 
especially those that are orthogonal with respect to the inner product 



(1-1) (/,<?>„ :=&m / f(x)g(x)W fl (x)dx, 

J« d 

where W^x) := (1 - \\xfY on B d , [i > -1, and & M is a normalizing constant such 
that (1, 1) = 1. The condition fj, > — 1 is necessary for the inner product to be 
well defined. The orthogonal polynomials for (/, g) are often called classical and 
they are eigenfunctions of a second order linear differential operator D^. 

The purpose of this paper is to study orthogonal polynomials on the unit ball 
with respect to an inner product that involves the integral of V/ ■ \7g. The study 
in this direction was initiated in [12], where the inner product 

(1.2) (f,g)j ■= — f V/(s) ■ Vg(x)dx + — f f(x)g(x)da, A > 0, 

and another one that has the last term in the right hand side of (/, g) 1 replaced by 
/(0)g(0) were considered. The orthogonal polynomials with respect to the inner 
product (/, g)j have the distinction that they are eigenfunctions of the second order 
partial differential operator D_±, which is the limiting case of when fi — ^ — 1 
([5]). For the inner product (/, g) I , the main term is the first term, the integral 
over M d , which however is zero if both / and g are constant, and the second term 
is necessary to make sure that (/, f)j = implies f — almost everywhere. 
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In the present paper, we first consider an extension of (/, g) j that has the weight 
function in the integral over M d and will construct an orthonormal basis for 
this extension. We then study orthogonal polynomials with respect to the inner 
product 
(1.3) 



f{x)g(x)W fl {x)dx + X / V/(a;) -Vg{x)W^x)dx 



A > 0. 

In this case, both the first term and the second term are over the unit ball and it is 
no longer clear which one is the dominating term. Nevertheless, we shall construct 
a sequence of orthogonal polynomials with respect to (/, g) M B d explicitly, which 
depends on a two-parameter family of Sobolev polynomials of one variable. The 
latter ones can be expressed as a sum of Jacobi polynomials with the coefficients 
computed by a recursive relation. The norm of the orthogonal polynomials with 
respect to {f,g)^ K d can also be computed by a simple recursive relation. Such an 
orthogonal basis can be useful in the spectral method for PDE [7] . 

In contrast to the case of one variable, the Sobolev orthogonal polynomials in 
several variables have been studied only in a few cases and mostly on the unit ball 
O [6l [9j [TTJ [12] . The setting of the unit ball has turned out to be more accessible 
and has revealed several interesting phenomena, which has stimulated and provided 
hint for the structure of the Sobolev orthogonal polynomials on other domains (see 
[2]). The present work is a further study on the unit ball, the results illustrate the 
impact of the secondary term in the Sobolev inner product, and also reveal further 
properties of the classical orthogonal polynomials on the unit ball. 

The paper is organized as follows. In the next section, we state the background 
materials and prove several properties of orthogonal polynomials on the unit ball 
and spherical harmonics that involve the action of V. Orthogonal polynomials 
with respect to the Sobolev inner product that extends (II. 2[) to weighted setting 
are discussed in Section 3, while those that are orthogonal with respect to (•,•}„ B d 
in (|1.3j) are discussed in Section 5. The Sobolev orthogonal polynomials of one 
variable that will be needed for the construction in the Section 5 are studied in 
Section 4. 



2. Orthogonal polynomials on the ball and spherical harmonics 

In this section we describe background materials on orthogonal polynomials and 
spherical harmonics that we shall need. The first subsection recalls the basic results 
on the classical orthogonal polynomials on the unit ball. The second subsection 
contains several lemmas on the spherical harmonics under the action of the gradient 
operator. The third subsection collects properties on the Jacobi polynomials that 
we shall need later. 

2.1. Orthogonal polynomials on the unit ball. For x, y £ M. d , we use the usual 
notation of and (x, y) to denote the Euclidean norm of x and the dot product 
of x, y. The unit ball and the unit sphere in R d are denoted, respectively, by 

M d := {x e M. d : \\x\\ < 1} and S^ 1 := {£ e R d : ||f|| = 1}. 

For [i £ R, let W M be the weight function defined by 

w„(x) = (i-\\xfr, \\ x \\<i. 
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The function is integrable on the unit ball if /j, > — 1, for which we denote the 
normalization constant of by 

h ( f w < \ri V 1 ^(/i + d/2 + 1) 

The classical orthogonal polynomials on the unit ball are orthogonal with respect 
to the inner product 

(2-1) (f,g) =b J f(x)g(x)W ll (x)dx, 

which is normalized so that (1, 1) = 1. 

Let H d denote the space of polynomials in d real variables. For n = 0,1,2,..., 
let denote the linear space of polynomials in several variables of (total) degree 
at most n and let V% denote the space of homogeneous polynomials of degree n. It 
is well known that 

^Iii = ( n + d \ and ^Vt=( n + d - l \:=rt 



A polynomial P G LT^ is called orthogonal with respect to on the ball if 
(P,Q)^ = for all Q G n^_ 1; that is, if it is orthogonal to all polynomials of 
less degree. Let V^iW^) denote the space of orthogonal polynomials of total degree 
n with respect to W^. Then dim V^iW^) = r%. 

For n > 0, let {P"(x) : \a\ = n} denote a basis of V^W^). Then every element 
of VniW^) is orthogonal to polynomials of less degree. If the elements of the basis 
are also orthogonal to each other, that is, (P™, Pg )„ = whenever a ^ (3, we call 
the basis mutually orthogonal. If, in addition, (P„,Pq) p = 1, we call the basis 
orthonormal. 

In spherical-polar coordinates x — r£, r > and £ G a mutually orthogonal 

basis of V^(W M ) can be given in terms of the Jacobi polynomials and spherical 
harmonics. 

Harmonic polynomials of d-variables are homogeneous polynomials in that 
satisfy the Laplace equation AY = 0, where A = + . . . + #3- is the usual 

1 d 

Laplace operator. Let denotes the space of harmonic polynomials of degree n. 
It is well know that 

Spherical harmonics are the restriction of harmonic polynomials on the unit sphere. 
If Y € then Y(x) = r n Y(£) in spherical-polar coordinates x — r£, so that Y 
is uniquely determined by its restriction on the sphere. We shall also use to 
denote the space of spherical harmonics of degree n. Let da denote the surface 
measure and uid denote the surface area, 

The spherical harmonics of different degrees are orthogonal with respect to the 
inner product 

</,0> S <-i:=— / /(0ff(0M0- 
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Let Pn (t) denotes the usual Jacobi polynomial of degree n, which is orthog- 
onal with respect to the weight function 

w a ,f}{t) = (l-t) a (l + tf, a,/3>-l, te[-l,l]. 

We can now state the orthogonal basis of V^(W^) in spherical-polar coordinates 
(see, for example, [5]). 

Lemma 2.1. For n <E N and < j < n/2, let {Y^ 2j : 1 < v < a d n _ 2j } denote an 
orthonormal basis for Hn-2j- Defir 



me 



(2.2) P*Jx) = P ^-V+^\ 2 \\xf - 1) Yr 2j {x). 

Then the set {P™ [/ (x) : 1 < j < n/2, 1 < v < af l _ 2 j} is a mutually orthogonal basis 
o/V„(W]u). More precisely, 

(P^x^P^x))^ = H£j n , m 5 jtk 6 U , V , 
where is given by 

(9 ~, rr^ _ (M+l)j(|)n-j(»-j+M+^) 

1 J J> ' il(M + ^)n-i(n + M+|) ' 

where (a) n = a(a + 1) . . . (a + n — 1) denotes the Pochhammer symbol. 

It is known that orthogonal polynomials with respect to W 7 ^ are eigenfunctions 
of a second order differential operator 2?^. More precisely, we have 

(2.4) XVP=-(n + cO(n + 2 M )P, VPeV> M ), 

where 



j=i j 



2/i + z t 



2.2. Spherical harmonics. For later study, we will need several properties of 
spherical harmonics. Let V denote the gradient operator 

V/:= (d 1 ,f,d 2 f,...,d d f) T , 

where dk denotes the partial derivative in the fc-th variable. We denote the dot 
product of V/ and Vg by V/ • Vg. 

In spherical-polar coordinates x = r£, r > and £ G § d_1 , the differential 
operators V and A can be decomposed as follows (cf. [3]): 

(2.5) V =iv +£|-, 

r or 

(2 6) A-— + — -— + — A 

dr 2 r dr r 2 

The operators Vo and Ao are the spherical parts of the gradient and the Laplacian, 
respectively. The operator A is called the Laplace-Beltrami operator, which has 
spherical harmonics as its eigenfunctions. More precisely, it holds(cf. [3]) 

(2.7) A Y(0 = -n(n + d-2)Y(0, VY G U d n , £ G S^ 1 . 

In addition, harmonic polynomials satisfy the following properties: 

Lemma 2.2. Let {Y™ : 1 < v < a„} be an orthonormal basis ofH^. Let x = r£, 
with r > and £ G § d_1 . Then 



SOBOLEV ORTHOGONAL POLYNOMIALS ON THE UNIT BALL 



5 



(i) £ • V Y?(x) = 0, 

i 71 Til 

(ii) WY-(x) ■ VY™(x) = ^ V Y?(x) ■ V Y v m (x) + —Y^{x)Y^{x). 

(iii) For 1 < v < and 1 < rj < a^, the following relation holds 

(2.8) — / VY?(0 ■ VY™(0da(0 = n(2n + d- 2)6 n , m 6„ >n . 

Proof. Since FJ 1 is homogeneous, 4^Y™(x) = nY™(x)/r and, by the Euler's equation 
for homogeneous polynomials, x ■ VY™(x) = nY™(x). Hence, by (12.51) 

77 

£ • V Y?(x) = x ■ VYl\x) - r-Y?(x) = nY?(x) - nY?{x) = 0. 

Using (i), the proof of (ii) follows from (I2.5[) and a straightforward computation. 
For the proof of (iii), we will need the Green's formula on the sphere (cf. [4]), 

(2-9) / Vo/(0 • V g(0da(0 = - [ A Q f(0g(0da(0. 

Applying (|2 .9[) on the identity in (ii) and using (|2.7I) . we obtain 

vrj l (e)-vr„ m (C)da(0 

VoC(0 • VoY^(OMO + nm [ Y, n (x)Y v m (x)da^) 

s-i-1 



= n(n + d-2) Y?{0Y™(0da(0 + n^ d 5 n , m 5^ 

= n(2n + d- 2)u)d5 n , m 5 1/tn . 
This completes the proof. □ 

Corollary 2.3. Let {Y™ : 1 < v < a^} be an orthonormal basis of H„. For 
/i > -1, n ^ m, 1 < v < and 1 < rj < a^, 

(i) b, [ Y?(x)Y™(x)W,(x)dx = ( f" 6 n , m 6v,r, 

JK d + t + 2 Jn 

(ii) b^ f vy;(!) • vY™{x)w^ +1 {x)dx = 2n(/x + 1)- — ( f )n . W^- 

Jn d 1 (/J+ 1 + 2)n 

Proof. The proof uses the following well known identity 

(2.10) / f(x)dx=f 1 r d - 1 f f(r0M0dr 

Jm d Jo Js d - 1 

that arises from the spherical-polar coordinates x = r £, £ e § d_1 . Since F„™ is a 
homogeneous polynomial, Y£{x) = r n Y^{^), so that 

b^ f Y^(x)Y^(x)W fl (x)dx = b M u d f 1 r 2n+d -\l - r^drS^S^ 
Jm d Jo 

by the orthonormality of Y™ with respect to (/, g) sd ~i, from which (i) follows upon 

evaluating the integral and using the fact that 



(2.11) b^ 



, £(m+i+|) 
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Since VY£(x) = ^'(Vy;)^), (ii) follows similarly. □ 

2.3. Classical Jacobi polynomials. We collect the properties of the classical 
Jacobi polynomials P^^\t) that we shall need below, most of which can be found 
in [TJ chapt. 22] and [TU]. The Jacobi polynomial pff ' (t) is normalized by 

(2.12) Pt ^ 1)= f- + -) = ^±lk > 

\ n J n\ 

For a, /3 > — 1, these polynomials are orthogonal with respect to the Jacobi weight 
function 

w a ^(t) :=(l-t) a (l + tf, -1<x<1, 
and their L 2 norms in L {w a g, [—1, 1]) are given by 
(2.13) 

- C p^Ht? w 3 (t)dt= 2 ° +P+1 r(» + a + i)r(» + /3 + i) 

n J_ x n [)WaA) 2n + a + (3 + l n\ T(n + a + + 1) ' 

The polynomial Pn'^\t) is of degree n and its leading coefficient k^ 1 '^ is given by 

(2.14) k^ = ±/ 2n + a + ^ 
The derivative of a Jacobi polynomial is again a Jacobi polynomial, 

(2.15) |p^)w = n+Q + /3+1 ei w) (o. 

The following relations between different families of the Jacobi polynomials can 
be found in (1, Chapt. 22]): 

(2.16) (1 + t)Pf +1) (t) = - gtL ^(t) + - J^— P^P(Jt), 

2n + a + p + 2 2n + a + p + 2 

(2.17) P^>(i) = a^P^\t) - 6 (f«.«p(°+ 1 .« ( t ), 

(2.18) P(^>(i) = a^P^+^t) + e-)p^f +1) (t), 
where 

f2 191 - n + a + f ) + 1 b <f*,p) _ n + P 

1 J} n 2n + a + p + V n ~2n + a + /3 + r 

They are used to prove the following two relations that we need later. 
Lemma 2.4. For a > — 1 and (3 > 0, 

(2.20) (1 + i) ^^(t) ^P^W+nP^ 1 ^ 1 )^) 

(2.21) ^ Q '«(*) + (1 + *) |^ a,/3) W = (0 + n) P^ 1 '^ 1 ^) 

Proof. From (|2 . 15[) and (|2.16l) . we deduce 

3 



(1 + t) = (1 + 1) n + a ^ + 1 Pi"-l W) ft) 

n + a + (3 + 1 



2n + a + (3 + 1 
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replacing the last term by, according to f|2 . 1 8[) . 

n + a + p + 1 n + a + p + 1 

and simplifying the result, we obtain (|2 . 20|) . The relation (I2.2ip is deduced from 
(I2~20l) upon using P^ l) ~ 1] {t) - P^~ 1J) \t) = P^f \t), which follows from (l2~T7l) 
and (|2~T8|) . □ 

3. SOBOLEV ORTHOGONAL POLYNOMIALS WITH A SPHERICAL TERM 

The first study of such orthogonal polynomials was initiated in |12) . where or- 
thogonal polynomials with respect to the inner product 

(3.1) </,<7) v := — / V/(x) • Vg{x)dx + — [ f(Og(OMO 

were studied. Let V„(V) denote the space of orthogonal polynomials of degree n 
with respect to (•, -) v . An orthogonal basis for this inner product is given as follows. 

Theorem 3.1. A mutually orthogonal basis {U™ v :0<j<t|,1<z/< a^_ 2 j} for 
V„(V) is given by 

KJx)=Y?(x), 

(3 2) 

u?a*) = (! - M 2 )Pj-r 2j+ ^\nx\\ 2 ~ wr 2j (x), i <j < % 

where {Y™- 2 3 : 1 < v < a*_ 2j } is an orthonormal basis of T-L^_ 2 j ■ Furthermore, 

(3.3) =nX + l, (U^, E/£„) v = A. 

In particular, it was observed that the space V^(V) can be decomposed in terms 
of the space V n (W^) of classical orthogonal polynomials on the ball and the space 
7i d n of the spherical harmonics; that is, for n > 1, 

(3.4) V*(V) = K © (1 - \\x\\ 2 )Vt 2 (Wi). 

Moreover, it was shown in [9] that the elements of V„(V) are eigenfunctions of the 
differential operator that is, they satisfy the equation (|2.4|) with /i — — 1. The 
phenomenon of the decomposition (|3.4j) was first observed in [11) for orthogonal 
polynomials with respect to an inner product that involves second order partial 
derivatives. 

The proof that U? are orthogonal with respect to (•, -} v given in [12] is based 
on the Green's identity 

Vf(x)-Vg(x)dx = f f{£)jLg(£) fa®- I f(x)Ag(x)dx. 
J§d-i dr J vd 

There is a weighted extension of Green's formula, which states that 

(3.5) [ Vf(x) ■ Vg(x) h(x)dx = [ /(£) ^ (0 fc(0 da 

.Ind ./sd-i dr 

f{x) (Ag(x) h(x) + Vg(x) Vh(x)) dx, 



s 
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where /, g, h are diffcrcntiable functions for which the integrals are finite. The proof 
of this identity follows from applying the divergence theorem on the identity 

div {fVgh) = Vf -Vgh + fAgh + fVg-Vh, 

where div is the divergence operator defined by div(/i(cc), f2(x), ■ ■ ■ , fd(x)) T — 
Ei=i<9i/i(a:). 

As an extension of (|3.ip . one can consider more generally the weighted inner 
product defined as follows. 

Definition 3.2. For /i > — 1 and f,g £ H d , we define 

(3-6) (f,g) :=—[ Vf(x)-Vg(x)W^ +1 (x)dx + —[ /(Os(0<MO- 

It is easy to see that this defines an inner product. The inner product (/, g) v in 
(|3.I[) is the limiting case when fj, = — 1. 

Let us denote by V^(V, W M ) the space of orthogonal polynomials of degree ex- 
actly n with respect to this inner product. We want to construct a basis for this 
space. 

Recall that the polynomials P™ l/ (x), defined in (|2.2|) . are orthogonal with respect 
to Wfj, on M d . To construct an orthogonal basis for V^(V, W^), we need the fol- 
lowing proposition that is of independent interest, which states that the orthogonal 
polynomials P^ v {x\ W^) satisfy another orthogonality on the ball. 

Proposition 3.3. Let fi > —I and let Pj l v (x) be the mutually orthogonal polyno- 
mials in V d {W^), defined in (|2.2[) . Then 



(3.7) 6„ / VPjUx) • VP^Jx^+^dx = W n {V)8 n , m 8j, k 8 v „, 

where, with HJ defined in (|2.3p . 

H^ n (V) := (n(2j + j(2j - d + 2)) H* n . 

Proof. For fi > —I, W / M+ i(x) = (f — ||a;|| 2 )W / A1 (x) vanishes when x 6 § d_1 . Hence, 
by (1331), 



(3.8) f VP^{x)-VP^(x)W^ +1 {x)dx 

P^(x) [(I - ||x|| 2 )A - 2(/x + l)(x, V)] P^W^*)^. 



Observe that [(I - ||x|| 2 )A - 2(/x + V)] has the same degree as Q, this 
proves the orthogonality when n^mby the orthogonality of P™^ e V^W^,). To 
prove the mutual orthogonality in (|3.7[) . we need to compute 

[(I-||x|| 2 )A-2( / i+f)(x,V)]P^(x). 

Let fa = n - 2j + Then P£ v (x) = P^ ] '(2r 2 - l)r n - 2j Y^- 2 ^ (£) in the 

spherical-polar coordinates x = r£. Using the expression of A in (|2.6[) and the 
derivative formula of the Jacobi polynomials, it is not difficult to verify that 

(1 - \\x\\ 2 )AP^(x) = 2(j + f , + l3 J + l) \2(j + » + ft + 2)r 2 P^+ 2A+2) (2r 2 I) 



-(2n - 4j + f f)P^+ 1 ' ft+1) (2r 2 - 1)1 Y^{x) 
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Furthermore, by $Tfy . x-V = £-V +r^; and, by (i) of LemmadJlJ £-V F™(£) = 0. 
It follows that 

(x, V)P"» = \2{j + fi + ftj + l)r 2 P^ 1 ' ft+1) (2^ 2 - 1) 



-(n - 2j)P. 



(eft), 



2r 2 - r 



Both of these terms are of the form gdlxlj 2 )!^™ 2j (x). Hence, looking at the leading 
coefficients of the Jacobi polynomials, given in (|2.14p . we conclude that 

[(l-\\xf)A-2( f ,+ l){x,V)]P^(x) 

= (n(2j + M + 1) - j(2j - d + 2)) P^ ] {2\\xf - l)Y?-*(x) + g(x), 

where g is a polynomial of degree at most n — 2. Consequently, by the orthogonality 
of P™, as an element of V^(W M ), (03) follows from (|5^| . □ 

When j = 0, the orthogonal polynomials Pj l v (x) become spherical harmonics 
Pqu( x ) — Y„{x), an d l|3.7p agrees with (ii) in Lemma l2~3l 



Theorem 3.4. For j < n/2, let {Y™ 2j : < v < a^_2j} be an orthonormal basis 
"fl ('!, >,■ Define 

(3 - 9) QUx):= [P^ ) (2||x|| 2 -1)-P^ ) (l)]yr 2 ^(a;), l<3<\- 

where /3j := n - 2j + Then {Q™„ : < j < n/2, 1 < v < a d n _ 2j } forms a 

mutually orthogonal basis of V^(V , W^) . Furthermore, 



r(/i + 2)r(n + |) 



r(/x + n + i + |) 

<Qi,„, Q?,„ W = A (»(2J + M + 1) - j(2j -d + 2)) 



+ j + l)r(n - j + f )(m + n - j + I) 



+ A(n - 2j) 



i!r( M + n-i + l + f)(n + M+|) 

r( A1 + 2)r(7i-2j + f)( M + i) 2 



r(n-2 J + M + f)r 



Proof. If j = 0, then the orthogonality of Qq v to lower degree polynomials and 
to Q™ v for j > follows from (i) in Lemma 12.31 and the orthogonality of Y„. 
Furthermore, the same consideration shows that 



LU d b 



— 2n(/x+l) — 



(/* + 1 + |)„ 



1, 



which simplifies to the stated formula upon using (|2.11[) . 

If < j < n/2, then Q™ u \§d-i vanishes on the sphere as its radial part vanishes 
when 1 1 a; 1 1 = 1, so that we only needs to consider the integral over M d in the inner 
product {•, -) v w . Furthermore, in terms of Pj\, in (|2.2[) of V d {W^), we can write 



Ux) = PUx)-P\ 



'{i)Yr 23 {x) 
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and observe that FJ 1 " 2 - 7 is a polynomial of degree less than n for j > 0, the mutual 
orthogonality of {Q^ v : < j < n/2} then follows from ()3.7|) . Moreover, it follows 
readily that 



+ ^2^ / [Yr 2j (x)?W» +1 (x)dx 



J- 

where we have used p^' n 2j+ 2 ■'(l) — i^ilii. Hence, the norm of u can now 
be deduced from (|3.7|) and (ii) of Lemma 12. 3[ and simplified to the stated formula 
upon using and (pETTj) . □ 

If we set /i — > — 1 in the above theorem, then since jTOl (4.22.2)] 

pi- M) w = ^(i-^lfw, 



it is easy to see that Theorem 13.41 agrees with Theorem 13.11 The decomposition 
(13.4[) . however, does not extend to V^(V, W M ). 

Finally, it is worth to mention that the orthogonal polynomials in CV^(W^) 
are automatically orthogonal with respect to a Sobolev inner product according to 
Proposition 13.31 



Theorem 3.5. Let \x > — 1 and let P™ v {x) be the mutually orthogonal polynomials 
in Vn(Wfj,), defined in (|2.2p . Then they are also mutually orthogonal with respect 
to the Sobolev inner product 

(3.10) [/,$]„:= ftp f f{x)g{x)W ll {x)dx + \ j V f(x) ■ Vg(x)W fl+1 (x)dx 

where A > is a fixed constant. 

Notice that the parameter of the weight function in the second integral is fi+1. 
It is possible to add terms with higher order derivatives with matching weight 
functions in the inner product, for which P™ v remain to be orthogonal. As we shall 
see in the Section 5, the orthogonal structure becomes far more complicated if we 
want the weight functions in the two terms of (|3.10p to be the same. 



4. A family of Sobolev orthogonal polynomials of one variable 

In the next section, we will construct an orthogonal basis for the Sovolev inner 
product (|1.3|) . As in the cases of (|2.2j) and (|3.9[) . the basis can be constructed in 
spherical-polar coordinates, where the main terms are of the form 

gi (2||x|| 2 -l)Y„"- 2 '(z). 

The polynomial qj of one variable, however, is rather involved in this case, which 
we now define. 
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Definition 4.1. Let A > 0. For a > — 1 and fj > max{0, ^y^}, we define a Sobolev 
inner product by 

(f,9)a,p ■= J f{t)g(t)w a ,p(t)dt 

(AY) ^0\[ 1 <ff'\( B(f3,d)(l+t)\ (g\ 

(4.1) + 2A/ (/,/)( md)(1 + t) 4(1 + i) a \l gl )w a , fi -i(t)dt, 



-l 



where A(/3, d) = /3B(/3, d) and B(J3, d) = 20 — (d — 2). 

In the case of j3 = ^=^, we have A(/3, d) = B(f3, d) = and the inner product 
takes a more familiar form 

(/,<?W = J J(t)g(t)w a ^(t)dt + 8X J f'(t)g'(t)w aif ) +1 (t)dt, 

that is, a Sobolev inner product associated to a coherent pair of measures (see [H]). 

We assume that /3 > max{0, ^j^} so that the 2x2 matrix in (|4.1I) is positive 
definite, which guarantees that {f,g) a ,/3 is indeed an inner product. Consequently, 
the orthogonal polynomials with respect to this inner product exist and they are 
uniquely determined up to a multiplicative constant. 

We let qj (t) denote the orthogonal polynomial of degree j with respect to 
{f,g) a ,p and normalize it so that it has the same leading coefficient as the Jacobi 
polynomial Pj a '^\t), that is, 

(4.2) q^ (t) = kf^ ] t j + lower degree terms, .^]_f 2 J + a + l 3 

With this normalization we have ?o ^(t) = Pq*' (t) = 1. 

Recall that the norm of the Jacobi polynomial p^ a '^ is denoted by hj^\ given 
in (|2.13j) . Let us denote the norm of qf"''^ in (•, -) a ,p by hj° , that is, 

j>>0) ._ ( '(.<*,?) J a >P)\ 

■ {''■> ■''■> ),,■ 

Our first result is to relate qj '^ with the Jacobi polynomials. 
Proposition 4.2. Let fi > -1 and f3> max{0, (d - 2)/2}. Then, for j > 0, 

(4.3) Pj»&(t) = a<f^q^\t) + d^\\)q^\t), 
where d^ a '^ are defined by 

(4.4) d^(A)=0, d^\x) = -b^ i ) J ( ^ Ty i=0,l,2,..., 



and aff 1 '^ and b^ 1 '^ are given by (|2.19|) . 



3 3 

Proof. Since {q\^ +1 '^\t)}k>o is a sequence of orthogonal polynomials with respect 
to (•, •) jU _)_i !/ g, we can expand P^^\i) in terms of them, that is, 

P ^\ t) ^W +m (t) with di = {P > ■ 

k=0 h k 
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Since q\^'^ is normalized so that its leading coefficient is the same as P^\ it 
follows from (1217)1 that = a ( f' 0) . 

We now calculate d J k for < k < j. From (|4.1j) . we need to compute 



P, 
-l 

2A 



A(/8,d) 



P, 



The first term in the right hand side vanishes for k < j — 1 since, by (|2.17l) and the 



orthogonality of P^ 



i r 



>>/8) pO-H-1,/3) 



«: / , ;"" ri, " > (t)-6y*"''pj-_ 1 



p(M+l,/3), 



(t) g^ +1 ^(*K +M Wdt = ! 



whereas for fc = j — 1, using (^"'-f (t) = Pj"f\t) + ... in addition shows that 

= -bf* f 1 pj^Ht)^f' p \t) w , +1 At)dt = -bf^hf_f^. 



An analogous reasoning, together with (|2.15|) . shows that the fourth and the fifth 
terms vanish for < A; < j — 1. Thus, we only need to compute the second and the 
third terms. This is where Lemma 12.41 is needed. Indeed, assume j3 > first; then, 
since A(/3,d) = /3B(/3,d), the relation (j2~2l"j) shows that, for < k < j - 1, 

A(f3,d) [ P^\t) q i^ ) (t)w tl+1 ^ 1 (t)dt 



■B(/9,d) 



1 d 



P 

_i<ft J 



(^,/3)/.n (M+1,/3) 



[2/?-(d-2)](/3 + j) 



P (^-i)( t ) ? ^i.« (t)u , /t+1 ^_ l(t)dt = 0j 



whereas if /3 = then A(/3, d) = B((3, d) = and these terms are automatically 



zero. 



Therefore, only the first term is nonzero and it is nonzero only for k = j — 1. 
Consequently, 

P\^\t) = a\^q ( f +m (t) + d^(X) q ^ U3 \t), 
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where djl'f^A) is given in (|4.4[) . This completes the proof. □ 

If we knew how to compute (^'^'(A), we could deduce recursively by 

(|4.3[) . However, since d) depends on the norm of q^ +h0) , this has a ring of a 

futile loop. It turns out, as our next proposition shows, that d^'^' can be deduced 
from a recursive relation. 

Proposition 4.3. The coefficients d^'f^(A) satisfy a recurrence relation 

(ak-\ H <Jt ' p) (X\ = 13 + 1 

{ ' ° {) (M + /3 + 3)[l + A(/i + /3 + 2)(/3-(d-2)/2)]' 

(4.6) d ( j fl,0) (X) = —- j-= —- , j>l, 

[BW + XC^ ) ] + df^ ) {X) 

where 

(4 7) a^ J3) = v±e±£± IK? + a + l M + & + iKgj + a + g) 

l ' J j j(2j + M + /3+l)(2j + M + /3 + 2)(2j + M + /3 + 3)' 

(4 - 8) - 1+ . ? (2.+M + /3 + 2)' 

(4.9) C?W> = (2 ' ? + i± + gMOj + l)(2/8 - (g - 2)) + 4j(j + /x + g + 1)] _ 
J 2j 

Proof. Since q^ +1 ^ 3 \x) = 1, the explicit expression for dg M '^(A) in f)4.5[) is deduced 
directly from flO}, upon using (|2"TT5|) . ([235]) . and the fact that 

h$j*W = = + 2A J 4(/?,d)^ +1 ^- 1 \ 

which follows directly from the definition of (•, •) M +i,^. 



For j > 1, we use (|4~3|) and (144)) to obtain 



9 f*,U u ~'~i,PA2 



\ 3 J 3+1 d (M,P)(x) 3 3 3 



Solving ei^'' 3 '' (A) from the above identity leads to 



df'«(A) = - V ; ' ' + l 



Jm.M b (, P) MP) 



which is the recursive relation (14.61) with 



Now, using the formula for hj +1 '^ in (]2 . 13[) and the expressions for and 
j n (|2.19[) . it is easy to verify that has the explicit expression in (|4.7|) . 
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In order to evaluate B { f' P) and ', we need to calculate (P^\ P^ P) )^+i,p. 

From (|4T1) . 



(4.io) (^W^)^ = y 1 («)) w^w* 



+ 2A 



1 

1 



2i?(/9,d) J P^\t)-P^\t) Wlx+1 , {t)dt 
4 £ (i p ? ,P) ®) ^ + x,0+i(t)dt . 



Since A is a parameter, the first term in the right hand side is the numerator of 
Bj^'^, whereas the expression in the square bracket is the numerator of Cj . 
Using (I2.17[) and the orthogonality of the Jacobi polynomials, we see that 



so that, by (pTTSl) and (|2TT^|) again, 

Cj + M + j8 + l)(? + H + l)(2j + m + 0) , J +P 



j(2j + ti + f3+ l)(2j + n + + 2) 2j + M + + 1 
j(2j + M + /3 + 2) ' 

For cj' 1 '' 3 '', we compute the three terms in the square bracket of (|4.10p . First we 
combine the first term and half of the second term, and use the relation (|2.21[) to 
deduce 



C x :=A((3,d) j ^ (p^\tj) Wll+1 ^ x {t)dt 

+ B(f3,d) £ P^\t)jP^\t) Wll+1< p(t)dt 



3 

J -I 

ju(M,/3) 

: W> W + 3) k yl 1 , p - 1) h? +1 ' p - 1) = W ~{d- 2W + J>? +M " 13 - 
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Next, by (|2.15[) and the orthogonality of the Jacobi polynomials, half of the second 
term becomes 

C 2 :=B(P,d) j_P^\t)jP^\t)w il+x , l) {t)dt 



.„(/x+l,/9+l) 



Finally, using (|2.15p . the third terms becomes 

C 3 := 4 £ (jLpM)^ 2 w ^ +1 ( t )dt = (j + p + p + lfhf_\ 



Combining these terms we obtain 



3 l(a»,/3) i>+1,/9) 

which simplifies to the formula (|4.9I) . □ 

The recurrence relation @2>) shows that d^' 0) (A) can be expressed as a continu- 
ous fraction. Consequently, we can express d^'^\x) in terms of a rational function 
of A whose numerator and denominator are, respectively, the (j — l)th and jth 
elements of a sequence of orthogonal polynomials. 

Definition 4.4. For j 6 No, define the polynomials r^' 13 ^ by rjf'^(X) = 1 and 
(4.11) r^f(A) = (C^X + B^)r^\x) - A<1'f V£f } (A), 

where we assume r_{ (X) — and, /or j > 1, ylj A ' , ' 3 - ) j m^i/ 9 ) ana > are d e ji ne ^ 

in (|4.7[) . (I4.8[) . (|4.9[) . respectively, whereas for j = 0, 

_ (/i + l)(/3 + l)(/i + /3) ( ^) (m+1Km+^) 

' ( A i + /3 + 2)( M + /? + 3)' • (/x + j9 + 2) ' 

_ (M + l)(M + /?)(2/3-(rf-2)) 
G ° - ~ 2 ■ 

For fi > -1 and /3>(d - 2)/2, the coefficients A^' P) , s] M,/3) and cj^ 3 are all 

positive for j > 0. Consequently, by the Favard's theorem, the polynomials r^ 1 '^ 
are orthogonal with respect to a positive linear functional. Since these coefficients 
are explicitly known, one naturally asks if it is possible to identify these orthogonal 
polynomials, say with some classical orthogonal polynomials. The numerator of 
the coefficient C^'^ in (|4.9[) contains a quadratic polynomial in j that does not 
factor into product of linear factors for independent parameters of /3, fi and d. This 
shows that the orthogonal polynomials are not hypergeometric type in general. 
In some special cases, such as 3(fi + 1) = 4(/3 — (d — 2)) or (3 = (d — 2)/2, the 
quadratic term does factor and the corresponding orthogonal polynomials could be 
of hypergeometric type. For our study of the Sobolev polynomials on the ball in 
the next section, (3 = n — 2j + for integers n and j with < j < n/2, the 
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factorization happens only when n— 2j = 0. Since this is a sidetrack from our main 
purpose, we shall not pursue this direction any further. 

For our study, it is sufficient to remark that the three-term relation in (|4.f f [) 

offers an effective way to generate rj M '^ . 

Proposition 4.5. For j — 0, 1, 2, . . ., the coefficients d^'^(A) in (|4.4[) satisfy 
(4.12) ^\x) = -A^%J§-. 

Proof. The case of j = follows directly from (|4.5I) . which allows us to deduce, 
from r^CA) - 1 and r { f' f3) (A) = C^X + B^\ the formulas for A^'^, 
and Cq*'^\ Assume that ()4.12|) has been established for integers up to j — 1. Then, 
by (USD, 

_ ~j 



(nit*,?) , yr&'Ph ^WlVi l A ) 



rf^(A) 



+ XC^)r^\x) - a££V£«(A) ' 
which is (|4.12l) . By induction, this completes the proof. □ 

Corollary 4.6. Let fi + 1 > -1 and ^>max{0, (d - 2)/2}. Then, for j > 0, the 
Sobolev orthogonal polynomials in (14. 2|) satisfy 

i J n(^) 

a j r j ( A ) i=0 °i 
where a^' 13 ^ are given in (|2.19[) . r^' ^(A) are given in (I4.11[) . and D 1 ^ 1 are defined 
by Dfr'V = (fx + P)(jjl + + 1) and, /or J > 1, 

r4 U x n 0x,« _ 2 J Q- + M + g + l)(2j + M + m^Hj 1)! 

1 ' J j " ( Ai + l) J (/3 + l) J 

Furthermore, the Sobolev norm of q^ 1 "^ is given by 

,( M ,/3), (m+i,/3) jM,/3) m 

mk, / (^+i,« ( M +i,/3)\ _ gj+i n ] r 3 +i y*) 

where bf' P \ h^ +1 ' p) and A<f' P) are given in 1(235]) , ([2331) a " rf 02]), respectively. 
Proof. Using (I4.12[) . the relation (|4.3I) can be rewritten as 
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For j > 2, it is easy to check that A^/af' 0) = / , so that the above 

identity can be written as 

where tff ,P \t) = D^' P) r^' P) (X)q i j ' 1+m (t). Furthermore, setting j = 1 in 
shows that the above relation also holds for j — 1 (which is how the value of D^f'^ 
is determined). In particular, it follows that %' it) = D^ ,l3 \ Summing up the 
telescoping sequence shows that +1 (t) can be written as a sum in terms of the 
Jacobi polynomials, which proves (I4.13[) . 

Finally, (PTT5)) follows directly from P~lf and P~T2l . □ 



5. SOBOLEV INNER PRODUCT ON THE UNIT BALL 

Recall that our main task is to consider orthogonal polynomials with respect to 
the inner product (|1.3I) . which we restate below. 



Definition 5.1. Let A > and /i > — 1. For f,g£ H d , we define 



(/>ff)v,W u , B d : ~ 



f(x)g(x)W li (x)dx + X / Vf(x) ■Vg(x)W^(x)dx 



It is easy to see that (■, -) v w Md is indeed an inner product. Let V d (V, W^,M d ) 
denote the linear space of polynomials of total degree n that are orthogonal to all 
polynomial of lower degree with respect to this inner product. It follows then that 
dim V^(V,W M ,B d ) =r*. 

An orthogonal basis for V^(V, W M , M d ) is given in the following theorem. 

Theorem 5.2. Let A > 0. For < j < n/2, let fij := n-2j + ^ and let q^ ,0i) {t) 
be the k-th Sobolev orthogonal polynomial associated with the inner product (•, O/^.ft • 
Let {Y™~ 2 i : 1 < v < a d l _ 2 j} be an orthonormal basis ofH d l _ 2 j- Define 

(5.1) R-Jx) := qf' Pi) {2\\xf - l)Yr 2j {x). 

Then the set {i?™ v (x) : < j < n/2, < v < af l2 j} is a mutually orthogonal 
basis of V d (V,W^,M d ). Moreover, 



Proof. We need to calculate, 

(i?^ ^(a;), R kri {x))y ^ w B d =6p 



+ A f VRl^-VR^^W^dx 



We consider the two integrals in the right hand side separately. For the first integral, 
using (|2.10p and the fact that FJ 1-2 - 7 is homogenous and orthonormal with respect 
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t° (') ')§,d-i, we obtain, 



= b^ d [ qj {2r 2 - l)q k {2r 2 - l)r 2 ^ +1 (f - r 2 Y8 n ^,m-2kS^ 



2n-2j+fi,+d/2+l 



qj (t) q k (t)(l~tr,(l + t) n - 2] + ^dt 



Jn—2j,m- 



-2k8v,T], 



where we have omitted the superscript of q y p' L 2 ' for simplicity, which we 
shall adopt in the rest of this proof. 

For the second integral, we observe that, by the product rule of differentiation, 



W„) = ?;(2|W 2 - \)l Xi Y?- 2 1{x) + qj (2\\x\\ 2 - l)8iY?-*(x), 
which implies that 

VRl„{x) ■ VR^Jx) = qj (2 \\x\\ 2 1) q k (2 \\xf - 1) VY^~^{x) ■ VY™- 2k {x) 

+ 4(m - 2fc)< Z ;-(2||x|| 2 - l) gfe (2|| 2 ;|| 2 - l)Y?- 2j (x)Y™- 2k (x) 
+ 4(n - 2j)^(2||x|| 2 - l)q' k (2\\x\\ 2 l)Y v n ~ 2j {x)Y™~ 2k (x) 
+ W\\x\\ 2 q' 3 (2 \\x\\ 2 l) q ' k (2\\x\\ 2 - l)Yr 2j (x)Y™- 2k (x). 



We now integrate above expression term by term and apply (|2.10[) . For the first 
term, we obtain, using (|2.8|) . that 



qj (2\\x\\ 2 - l) qk (2\\x\\ 2 - l)VYr 2j {x) ■ VY™- 2k {x)W^{x)dx 



qj (2r 2 - l) qk (2r 2 - 1)(1 - r ^ r n-2j+m-2k+d-l 
(n-2j)(2(n-2j) + d-2) 

X o drbJdO n -2j,m-2kOv,n 



qj (t) qk {t){i-ty{i+t) n -v +d ' 2 - 2 dt 

i 

(n-2j)(2(n-2j) + d-2)u J d x 
X 2 n - 2 ^+^+ d / 2 °n-2j,m-2kOu,ri- 

The second and the third terms are similar. For the second term, we obtain 
4(m - 2k) f J, (2\\x\\ 2 - 1)%(2|M| 2 - 1)Y^ (x)Y^- 2k (x)W^(x)dx 



4(m-2fc) / ^(2r 2 - l)g fc (2r 2 - 1)(1 - r 2 )^r 
Jo 



drLOd8n-2j,m~2k5v,ri 



2 „-2j+ M +d/2 °n-2 J ',m-2fe5 I /,7 ? - 
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Finally, for the fourth term, we obtain 

16 / \\x\\%{2\\xf - 1K(2|M| 2 - {x)Y™-™{x)W^x)dx 

=16 f 1 q' J (2r 2 - l)q' k (2r 2 - 1)(1 - r 2 Yr 2n -^ +d+1 druj d 5 n ^ m . 2k 5 v , n 



m^m-tY(\+t) n - 2 i+ d / 2 dt— 



T $ n — 2j ,m~2k0~ v ,n - 



; | ' 2 n-2j+/t+d/2-2 

Putting all four terms together, we conclude that the second integral satisfies 

x (n - 2j)(2(n - 2j) + d - 2) | ^(*)« fc (t)(l - t)"(l + t) B_2, ' +<,/2 ~ 2 * 
+ 2(n - 2j) / (*)<&(*))'(! - *)" (1 + t)"-^/ 8 " 1 * 



Together, the first and the two integrals lead to 

(Kj,ui R T,n)v,W,,,B d = 2 n-2j+ij,+d/2 ^n-2j,m-2k^v, v {<lj i 5fc) /Jin _ 2j '+^2 ■ 

Using (12. lip , this completes the proof. 



□ 



By (|4.15p . the norm of can be expressed in terms of rf' , which can be 



computed recursively. Moreover, as a corollary of (|4.3[) . we also have the following 
relation between classical orthogonal polynomials Pj\, on the ball and our Sobolev 
orthogonal polynomials. Let us denote P™ v in V^W 7 ^) by Pf„{-\ W a ) and, similarly, 
denote RJ v (x) in V r f(V, W M ,B d ) by RJ v (x; W^) to emphasis the dependence on \i. 

Corollary 5.3. Let (j, > -1. For < j < n/2, let fa ■ = n - 2j + 4=2. . Then 

RfJx; W„) = a'f'^mjx; W M+ i) + (X)m: 2 Jx; W u+l ). 



Furthermore, 



2 J ( 



where the notations are same as those in (|4.13[) . 
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